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^ ■ Abstract 

> ; 

O . We show that two-photon transport can be modulated by a two-level emitter coupled to a 

^ I cavity in a one- dimensional waveguide. In the ordinary case, the transmitted light has a 

OO I wider frequency spectrum than the situation without the cavity because it is reflected and 

^ ' scattered many times. But when the two photons are resonant with the cavity resonance 

reflection frequency, the frequency spectrum of the transmitted light becomes narrower than 

that without the cavity. This means that properly tuning the cavity resonance frequency 

can improve the photon-photon interaction. In addition, we show that the two-photon 

C ■ intensity correlation functions are nearly opposite to each other at the two sides of the 

^ ' emitter transition frequency rather than be the same, which is exactly the Fano resonance 

O^. line shape for two photons. Such an effect is important for lowering the power threshold in 

optical bistable devices and for sensing applications. When the emitter transition frequency 

equals to the cavity resonance frequency for a high-Q cavity, our results agree with the 

■^ ' recent experiments and theories. 

m '. 

Q^ ■ Keywords: one-dimensional waveguide, atom-cavity system, two-photon scattering, Fano 

^ J ! resonance line shape 

!>■ 

o; 

1/ • 1. Introduction 

^ ' Quantum internet is important for quantum computation, communication and metrology 

j^ ■ [1]. The main methods to realize the quantum networks are the cavity- QED-based protocol 

[2] and the protocol provided by Duan, Lukin, Cirac, and ZoUer (DLCZ) [3, 4]. The former 
one uses the interaction between light and A-type three-level atoms in a high-Q cavity to 
store quantum states [5]. The latter one involves the measurement-induced entanglement 
[6], replaces single atoms by atomic ensembles so that it reduces the restriction of the 
cavity's quality and utilize the built-in entanglement-purification to create the long-distance 
quantum communication efficiently. 
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(a) Experimental setup (b) Two-level emitter 

Figure 1: (a) The schematic experimental setups for concurrence measurement of the correlated transmission 
probability density \ti{x\^X2)\^ . /?i, I?2 are photon-detectors with adjustable positions. BS is beam splitter. 
The "=" denotes the two-level emitter in the ID waveguide. The black parentheses symbolize the cavity, of 
which the length is 11. (b) The two- level emitter, k and p arc the momenta of the two photons. 51 is the 
emitter transition frequency. 

Thus it is of great significance to consider the interaction between few photons and 
a cavity-coupled emitter. Here we focus on a unique one-dimensional (ID) waveguide for 
photons because it is helpful to realize the strong photon-photon interaction [7, 8]. Recently, 
people have already done some researches about the interaction between light and an atom- 
cavity system in a ID waveguide both in experiment [9, 10] and in theory [11-17]. These 
theoretical approaches include either scattering matrix [11-14] or quantum field theory [15- 
17]. The two methods get the similar results for single photons [14, 15]. 

In this paper we study the interaction of two photons and a cavity-coupled two- level emit- 
ter in a ID waveguide. By analyzing two-photon wave packets, generalizing the scattering 
matrix [18] to dispose frequency-spectrum transformation, we get the two-photon correlated 
momenta distribution, correlated transmitted coefficient, Fano resonance line shape and 
other important properties. Comparing with most existing works, our results are more gen- 
eral because our approach is not restricted to a high-Q cavity. Because our model includes 
the decay inherently, our results exactly agree with the experimental results in Ref. [10]. 
In Sec. 2 we show the model and analyze two-photon overlaps. In Sec. 3 we generalize the 
scattering matrix [18] to deal with the frequency-spectrum transformation. We show the 
results in Sec. 4 and make a summary in Sec. 5. 

2. The Model and Two-Photon Wave Packets 

There are several methods to realize the atom-cavity system in a ID waveguide [9-11]. 
Fig. 1(a) shows the schematic experimental setups [8]. The incident beam is monochrome 
plane wave. After the two photons are scattered by the whole system, we separate them by 
a beam splitter (BS) and use two single-photon counters to detect them. We register the 
probability of concurrent detections. Similarly with the Hanbury Brown- Twiss effect [19], 
the results can show correlation of the two photons. Fig. 1(b) shows the parameters of the 
two-level emitter. 

Fig. 2(a) shows wave packets of the two photons. After the incidence of the two photons 
on the reference plane, they can be scattered to different sides of the plane. Incident photons 
can also come from different sides of the system. So we cannot divide these wave packets 
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Figure 2: (a) The ellipses are the two-photon wave packets, of which some are at together and some are 
separated. The arrows in the ellipses represent the moving directions of the photons. The uncircled numbers, 
e.g. "1" and "0" , denote the known amplitudes of wave packets each. The black, green and purple lines are 
the reference plane of the cavity, the free space and the emitter respectively, (b) Symbols in Eq. (1). Siur 
et al. are frequency spcctrums of the wave packets. The blue line is the reference plane. 

into left-coming or right-coming parts. However, we find that if each one of the two photons 
runs out of the cavity, it cannot form the situation that two photons are at together again. 
So we take this situation as "loss" at first [see the "uncoupled" wave packets in Fig. 2(a)]. 
After calculation of wave packet 0, we add wave packets (2) and (3) to make the complete 
resolution of this problem. 



3. Scattering Matrix and Frequency-Spectrum Transformation 

For two-photon case, the emitted light has a frequency-spectrum distribution [7, 8]. So 
we need to generalize the scattering matrix [18] to dispose this transformation. In fact, we 
need to solve several integral equations to get the spectrums of the outgoing wave packets 
[wave packets 0, and in Fig. 2(a)] from the incident wave packets [the wave packets 
denoted by uncircled numbers "1" and "0" in Fig. 2(a)], of which the amplitudes are known. 
Here we use scattering matrix to express and solve these equations. We define the scattering 
matrix as follow 



SiLL{k,p) 

S2RRik,p) 
SRLout{k,p) 



i?22 T22 T12 
T22 R22 T12 
T21 T21 Tii_ 



SlRR{k,p) 
S2LL{k,p) 
SRLin{k,p) 



(1) 



Here, SiRR{k,p) et al. are frequency spectrums of the wave packets [Fig. 2(b)]. The operator 
T22 (-R22) represents that two photons come together and transmit (reflect) together. T12 
describes that two photons come from different sides of the reference plane and are scattered 



to the same side, of which the reverse process is denoted by T21. Tn denotes that two photons 
come from and to different sides of the system. According to [8] we have 

1 f°° 

T22f{k,p) = - dkidpif{ki,pi)x 

{h^ipAKki - k)6ip, -p) + 6{k^ -p)6{p, - k)] + B6iE, - E)}, 

1 f^ 
R22f{k,p) = - dkidpif{ki,pi)x 

{fk,rp, m, - k)5{pi -p) + 5{h - p)Kpi - k)] + B5{E^ - E)}, 

1 . I r°° 

Ti2f{k,p) = —T2if{k,p) = - dkidpif{ki,pi)x 

[6{k^ - k)5{pi -p) + 5{k^ - p)5{pi - k)] + B5{Ei - E) 



''fci'^pi T f'ki'^pi 



(4) 



1 f°° 

Tuf{k,p) = - dkidpif{ki,pi)x 

ft_hh^l+^hIj!L[s^ki - k)Sipi -p) + Siki - p)Sipi - k)] + y/2BS{E, - E) 



(5) 



where f{k,p) = f{p,k) denotes the frequency spectrum of the incident photons, ik and f^ 
are the transmission and reflection amphtude of the emitter for single photons respectively 
[8]. B is the background fluorescence 

4zr2 Ei-2n + iT 

~ vr [4A2 - {El - 2fi + iVy] [AA^ - {Ei - 2n + iVy] ' ^^^ 

Here, F is the coupling strength between the photons and the emitter. 

From Eq. (2) we know that operator T22 has two parts: the one is multiplied by [6{ki — 
k)S{pi — p) + S{ki — p)S{pi — k)] and the other is multiplied by 6{Ei — E). The calculation 
of operators can be considered by each part respectively. In Tab. 1 we have the operator- 
multiplying formula (see Appendix A for details) 

Table 1: Operator- multiplying formula 



Operators 


5{Ei - E) 


6{ki -k)6{pi -p) 


T 


BoBiiAi)B2iA2) 


U{A) 


R 


CoCi(Ai)C2(A2) 


rr(A) 


RT 


BoBi{Ai)B2{A2)rr{A2) 
+ CoCi(Ai)C2(A2)tt(Ai) 

+ BoCoBi{Ai)C2{A2)f,nt 


tt{A)rr{A) 



where 



1 r+oo 

Unt = 2 / B2{y)Ci{y)dy. (7) 
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For calculation, we need the inverse operator, which satisfies 

f-'ff{k,p) = ff-'f{k,p) = f{k,p) = l[f{k,p) + f{p, 



1 



dkidpif{ki,pi)[6{ki - k)5{pi - p) + 5{ki - p)5{pi - k)]. 



Then we derive Tab. 2, which shows the operator-inversion formula (see Appendix B for 
details) 



Table 2: Operator-inversion formula 



Operators 



5{Ei - E) 



6{ki - k)6{pi -p) 



T 
f-i 



5o5i(Ai)52(A2) 

l+BpgM ttjAi) ttjA^) 



tt{A) 

1 

U(A) 



where 



Qint „ 



Biiy)B2iy) 
tt{y) 



dy. 



Having the formulae in Tab. 1 and 2, we can put the reduced scattering matrix 
emitter as follow (see Appendix C for details) 



SiLL{k,p) 

S2RR{k,p) 



R T 

f R 



SiRBik,p) 

S2LL{k,p) 



where 



r'2fi2(l-r'2fn)~'f2i. 



T = T22 

^ = i?22 + r'2f 12(1 - r'2fn)-'T2i. 
Here we make the substitution 



r" = v2e"'""'' X r' 



„'2 



^2iEl 



where r is the refiectivity of the cavity, 2/ is the length of the cavity. 

Similarly with the approach for single photons in Ref. [11], we rewrite Eq. (10) as 



S2RR{k,p) 
S2LL{k,p) 



f-RT'^R kf-^ 



SiRR{k,p) 

SlLL{k,p) 



(9) 
of the 

(10) 

(11) 

(12) 

follow 

(13) 



The scattering matrix of the cavity and the free space are 



and 



S2RRik,p) 
_S2LL{k,p) 

S2RR{k,p) 
S2LL{,k,p) 



4.2 r_ r2_ 

^ t2 i2 

J_ 

i2 
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SiRRik,p) 

SlLL{k,p) 

SiRRik,p) 

SlLL{k,p) 



(14) 



(15) 



We combine the scattering matrix of the cavity, the free space and the emitter in sequence 
to get the scattering matrix of the whole system 



T. 



J.Z r_ 722 

t' ,9 ±0 



42 



AEl 









-iEl 



T-RT-^R RT-^' 



AEl 









-iEl 



±Z £_ ]_ 






(16) 



Next, the procedure of our approach is described as follow [below the wave packets are 
shown in Fig. 2(a)] 

1. By using the formulae in Tab. 1 and Tab. 2, we derive the operators T~22 and Tsi2T~22- 
So we get the frequency spectrums of wave packets (T) and (2), which are T~22fin{k,p) 
and Ts^i2T~22fin{k,p) respectively [11]. Here fin{k,p) is the frequency spectrum of the 
incident two photons, which is a delta function and denoted by the uncircled number 
"1" in Fig. 2(a). 

2. By using the scattering matrix of the cavity [Eq. (14)], we derive the frequency spec- 
trums of wave packets (2) and ® from and wave packet @ from (|). 

3. By using the scattering matrix of the emitter [Eq. (1)], we derive the frequency spec- 
trum of wave packet (7) from (5) and (B). When one photon of wave packet (7) transmits 
the cavity, we get wave packet (3) from (7). 

4. we calculate the two photon transmitted wave packets associated with (2) and (3). In 
this situation, the problem turns to single photon transport, so we can use the results 
in Ref. [11] directly. 

5. By adding the two photon transmitted wave packets associated with (2) and (3) to 
the wave packet (T), we derive the frequency spectrum of the two photon transmitted 
wave packet finally. After fourier transformation, we get the two-photon wave function 
t2{xi,X2) in part of the out states, in which both photons are transmitted. 



4. Results 

4-1- Situation for Low- Q Cavity 

Here we consider a low-Q cavity. We set the cavity reflectivity r = 0.9, cavity transmis- 
sion amplitude t = i\/l — r^, the coupling strength of the emitter F = 0.004(27rc//). 

Firstly, We consider the situation that two photons are near resonant with the emitter 
transition frequency Q [i.e. 6E ~ 0, A ~ 0, see Fig. 1(b)]. We find that, in ordinary case 
[e.g. -E'i/2 = 0.240(27rc//), see Fig. 3(b)], the correlated transmission probability density 
1^2(3^1, 3^2) P contracts, which indicates that the transmitted light has a wider frequency 
spectrum than the situation without the cavity [7, 8]. This is because the photons reflected 
and scattered many times. But when the two photons are resonant with the cavity resonance 
reflection frequency [e.g. -E'i/2 = 0.125(27rc//), see Fig. 3(c)], the two-photon wave packet 
extends. So the cavity makes the frequency spectrum of the transmitted light narrower than 
that without the cavity. As a consequence, we achieve a strong photon-photon interaction 
by properly tuning the resonant frequency of the cavity, which is important for realization 
of controlled optical nonlinearity. 
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Figure 3: Two-photon correlated transmission probability density |i2(a;i,a;2)p. Two photons are near res- 
onant with the emitter transition frequency f2. (a) shows the situation without the cavity, (b) and (c) 
represent the influences of the cavity. The abscissa is the distance difference a; = xi — X2 of the two detectors 
(L>i, D2) [Fig. 1(a)], which is scaled by l/F. 



Then we study the non-resonant case that two photons are not resonant with the emitter 
transition frequency VL = 0.325(27rc)// (Fig. 4). Here we consider the intensity correlation 
function g^'^\xi,X2), which is defined as follow [19] 



9''Kx,,x,)=g^'Kx) 



(2)1 



{a+{xi)a+{x2)a{x2)a{xi)) 1^2(2:1, Xa)^ 



{a'^aY 



D 



:i7) 



where x = Xi — X2. -D is the normalization constant and independent of x. 

We show that the intensity correlation functions g'^'^\xi^X2) are nearly opposite to each 
other at the two sides of the emitter transition frequency Vt rather than be the same. 
The transmitted two photons exhibit bunching behaviors and super-Poissonian statistics 
{g^'^\Q) > 1) when E/2 < Q. (the upside figures of Fig. 4) and anti-bunching behaviors 
and sub-Poissonian statistics (^''■^^(O) < 1) when E/2 < Q (the underside figures of Fig. 4). 
We already know that, for single photons, Fano resonance exhibits a sharp asymmetric line 
shape with the transmission coefficients varying from to 100% over a very narrow frequency 
range [12]. So here we show the Fano resonance line shape for two photons. Because of its 
sharp asymmetric line shape, such an effect is important for lowering the power threshold 
in optical bistable devices and for sensing applications [11]. 

4 ■ 2. Situation for High- Q Cavity 

In a high-Q cavity, we can compare our results with the recent works done by Birnbaum 
et al. [10] and Shi et al. [17]. We set our parameters as follow to fit their conditions: 
the cavity's refiectivity r = 0.9996 and transmission amplitude t = i\/\ — r^, the emitter 
transition frequency fi = 2t\cJI (which is also the cavity resonance frequency when r is a real 
number), the coupling strength F = 0.002(27rc//). We also consider the intensity correlation 
function g^'^^{x\,X'i) = g^'^\x) (Fig. 5). 

Fig. 5(a) represents the relation between log[g^'^\0)] and E/2, which agrees with the 
results in Ref. [10, 17]. When E/2 = 0.9966(27rc//) or 1.034(27rc//), log[g^^\0)] gets its 
minimum. Comparing with Ref. [10, 17], we also consider the situation when E/2 = 
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Figure 4: Intensity correlation function g^'^'(x\^X2). Two photons arc non-resonant with the emitter tran- 
sition frequency il = 0.325(27rc)/L The incident two-photon momenta difference A = [see Fig. 1(b)]. The 
abscissa is the same as that in Fig. 3. 

0.9966(27rc//) [Fig. 5(b) and Fig. 5(c)]. Because our model includes the decay of the cavity 
inherently, our results exactly match the Figure. 3 of Ref. [10], except where has a modu- 
lation arising from center-of-mass motion of the trapped atom. In Ref. [17], however, the 
g^^^r) arises to unity after infinite time because it does not consider the decay of the cavity 
or that of the two-level emitter. 

5. Summary 

In conclusion, we study the two-photon correlated scattering by a two-level cavity- 
coupled emitter in a ID waveguide. Our results agree with the recent works [10, 17]. Our 
research is important for realization of strong photon-photon interaction at the case of few 
photons and for quantum network. In addition, we propose a general approach to two- 
photon, cavity-coupled and multi-frequency scattering problem in the ID waveguide. The 
same approach can deal with the three-level emitter, which is useful to make the photon- 
transistor [20, 21]. 



6. Acknowledgments 

We acknowledge the financial support of Special Prophase Project on the National Basic 
Research Program of China (No.2011CB311807), the National Natural Science Founda- 
tion of China (No. 10804092), the Natural Science Basic Research Plan in Shaanxi Province 

8 





0.5 



_JL^ 



0.992 0.996 1.000 1.004 1.008 
E/2 

i^)log{g(^\0))~E/2 



-1 1 



{h)g(^^{x)-x 



-20 -10 10 20 

X 

(c)5(2)(a;)-x 



Figure 5: Intensity correlation function 5^^^ (a;) for a high-Q cavity. The abscissa of (a) is scaled by {2ttc/1). 
The abscissae of (b) and (c) are the same as that in Fig. 3. Here (b) is different from (c) in the range of the 
abscissa. 



of China (No.2010JQ1004), the Fundamental Research Funds for the Central Universities 
(No.xjj20100099). 

Appendix A. Derivations of Operator-Multiplying Formula 

Here we derive the Operator- multiplying formula in Tab. 1. If operators T and R are 
defined as follow 



9{k2,P2) =ff{k,p) 



dkidpif{ki,pi)x 



{tt{k,,p,)[6{k, - k2)6{pi-p2) + 5{h-p2)5{p, - k^)] + fi(Ai, A2)5(^i - E2)}, 

1 f°° 
Rg{k2,P2) = - dk2dp2g{k2,P2)y< 



(A.l) 



(A.2) 



{rr{k2, P2) [5{k2 - k)6{p2 -p) + 6ik2 - p)5{p2 - k)] + C(A2, A)5(i?2 - E)}, 
then the operator RT can be derived as follow 

r"00 ^00 



RTf{k,p) 



^0 



dkidpidk2dp2f{ki,pi)x 



{tt{ki,pi)rr{k2,P2)[S{ki - k2)5{pi - P2) + 5{ki - p2)S{pi - k2)]x 

[S{k2 - k)S{p2 -p) + 5{k2 - p)5{p2 - k)] 

+ 5(Ai, A2)rr{k2, P2) [^(^2 - k)5{p2 - p) 

+ 5{k2-p)6{p2-k)]5{E,-E2) 

+ C{A2,A)tt{k^,p,)[5{k, - k2)6ip,-p2) 

+ 6ik,-p2)6ip,-k2mE2-E) 

+ 5(Ai, A2)C(A2, A)5{E, - E2)5{E2 - E)} 



2./0 



dkidpif{ki,pi)x 



{tt{ki,pi)rr{ki,pi) [6{ki - k)5{pi - p) + 5{ki - p)5{pi - 
+ [5(Ai,A)rr(A;,p) + C(Ai,A)rr(A;i,pi) 

/oo 
5(Ai, A2)C(A2, A)(5(E2 - E,)dA2]5{E, - E)} 
-oo 

dkidpif{ki,pi)y. 



oo 
oo 



(A.3) 



2.0 

{tt(A;i,pi)rr(/i;i,pi)[5(A;i - /i;)5(pi - p) + 5(A;i ~ p)5{pi - k)] 

+ [So5i(Ai)S2(A2)rr(A) + CoCi(Ai)C2(A2)tt(Ai) 

+ EoCo5i(Ai)C2(A2)/.nt]5(^i - ^)}. 

Here we factorize the background fluorescence i?(Ai, A2) and C(Ai, A2) as follow 

5(Ai,A2) = 5o5i(Ai)i?2(A2), 

C(Ai,A2) = CoCi(Ai)C2(A2). ^ ■ ^ 

For example, if T = T22 [Eq. (2)], then 

AiT^ E1-2VL + iV 



5(Ai,A2) 



TT [4A2 - [El - 21] + iVf] [4A2 - {El - 2fi + zr)2] 

4zr2 

5o = (^1 -2Vt + iV), 

IT 
52(A2) ^ 



In Eq. (A.3) 



4AI - {El - 2Q + iVy 



1 ^+00 
f^nt = - B2{y)Ci{y)dy. (A.5) 



In Eq. (A.l), we find that an operator can be determined by the part multiplied by 
[6{ki — k)6{pi—p) + 6{ki—p)6{pi — k)] (e.g. tt{ki,pi) in operator T) and the part multiplied 
by S{Ei — E2) (e.g. -B(Ai, A2) in operator T). So we briefly write the result of Eq. (A.3) in 
Tab. 1. 

Appendix B. Derivations of Operator-inversion Formula 

We use the operator-multiplying formula to derive the operator-inversion formula in 
Tab. 2. We order that R = T~^, so we have 

Rff{k,p) = fRf{k,p) = f{k,p) = hf{k,p) + f{p, k)] 

1 r ^ (B.l) 

= - / dkidpif{ki,pi)[6{ki - k)5{pi - p) + 5{ki - p)5{pi - k)]. 
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Comparing Eq. (B.l) with Eq. (A. 3), we get 

tt{ki,pi)rr{ki,pi) = l, 
5ofii(Ai)fi2(A2)rr(A) + CoCi(Ai)C2(A2)tt(Ai) + fioCoBi(Ai)C2(A2)/.„i = 0. ^ ' ^ 

So 



1 + Bof, 



int 



C.(AO . 1^. ,B.3, 



tt(A2) 

Here 



/.ni=2y B,{y)C,{y)dy = - j -Ll^±-^dy ^ g,^,. (B.4) 

We briefly write the result of Eq. (B.3) in Tab. 2. 

Appendix C. Derivations of the Reduced Scattering Matrix of the Emitter 

Bellow we derive Eq. (10) and Eq. (11). From the third equation of Eq. (1) we have 

SRLout{k,p) = f2lSiRR{k,p) + T2i5'2Ll(/c,p) + f2lSRLin{k,p). (C.l) 

By using the following relation (see Fig. 2(a)) 

SRLin{k,p) = r''^SRLout{,k,p), (C.2) 

we have 

SRUn{k,p) = r'\l - r'^Rur%i[SiRR{k,p) + S^LLik^p)]. (C.3) 

By putting Eq. (C.3) into the first and second equations of Eq. (1) 



SlLL{k,p) = R22SlRR{k,p) +f22S2LL{k,p) + f 12 SRLin{k,p) , 
S2RR{k,p) = f22SiRR{k,p) + R22S2LL{k,p) + fi2SRLin{k,p), 

we get 

r = f22 + r'2fi2(l-r'2i?i,)~'T2i, 

^ = i?22 + r'2f 12(1 - r'2i?n)-'T2i. 

Here we make the substitution 



(C.4) 



(C.5) 



r 



/2 



726^*^' X r^, (C.6) 



where r is the reflectivity of the cavity. 

By using the formulae in Tab. 1 and Tab. 2, we get the expressions of T and R. 
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